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Abstract 

We study r — > [iVifi, i = e,fi, r decays in the two Higgs doublet model, with the 
inclusion of one and two spatial non-universal extra dimensions. We observe that the 
branching ratio is sensitive to two extra dimensions in contrary to a single extra dimension. 
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1 Introduction 



Lepton flavor violating (LFV) interactions are rich to analyze since they exist in the loop level 
and they are sensitive the physics beyond the standard model (SM). On the other hand they are 
clean theoretically because they are free from the nonperturbative effects. The experimental 
work done and the numerical results obtained stimulate the theoretical studies on these decays. 
The fi — > e7 and r — > /ry processes are among the LFV decays and the experimental the current 
limits for their branching ratios (BR) have been predicted as 1.2 x 10~ n flj and 1.1 x 10 -6 [2|, 
respectively. 

There is an extensive theoretical work done on the LFV decays in the literature. Such 
interactions are studied in a model independent way in in the seesaw model ^j, in the 
framework of the two Higgs doublet model (2HDM) j2HH], in supersymmetric models [3 El El 

ma [mini uni- 

In the present work we study r — > // Pi % — e, //, r decay in the model version of the 2HDM 
with the inclusion of the non-universal extra dimensions. This process exists at least at one 
loop level in the model III. The lepton flavor violation is driven by the internal scalar bosons 
h° and A and the transition r — > is obtained. Furthermore the internal Z boson connects 
this transition to the output (see Fig. Q). Notice that we respect the assumption of the 
non-existence of Cabibbo-Kobayashi-Maskawa (CKM) type matrix in the leptonic sector and 
vanishing charged Flavor Changing (FC) interactions. 

With the inclusion of additional dimensions there may be an enhancement in the BR of the 
decay under consideration. The extra dimensions have been studied in the literature extensively 
[Ti]-|31j. These new dimensions could not be detected at present and the most favorable 
description is their compactification on a surface with small radii. From the 4D point of view 
this compactification results in Kaluza-Klein (KK) modes of the particles with masses regulated 
by the parameter R, which is a typical size of the extra dimension. If all the fields are accessible 
to the extra dimensions, the extra dimensional momentum, and therefore, the KK number at 
each vertex is conserved. Such extra dimensions are called universal extra dimensions (UED) 
and they are studied in various works [T3J [HE H3 H3 1201 I2J1 1221 I2H] in the literature. If the extra 
dimensions are accessible to some fields but not all in the theory, such type of extra dimensions 
are called non- universal extra dimensions (NUED) (see for example |23j 1261 I3UI I31j ). In this 
case the KK number at each vertex is not conserved and tree level interaction of KK modes 
with the ordinary particles can exist. 

In our work we take the extra dimensions as non-universal and assume that the first Higgs 
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doublet and the gauge fields are accessible to the extra dimensions, however the leptons, the 
quarks and the second Higgs doublet, which contains new Higgs particles, live in the 4D brane. 
In this case the contributions due to the additional dimensions comes from the internal Z boson 
KK modes after the compactification of the external dimensions on orbifold S 1 /Z 2 ((S 1 x 
S 1 )/Z 2 ) for a single (two) extra dimensions. We observe that there is almost two orders larger 
enhancement on the BR of the process we study in the case of two NUED. This enhancement 
is obviously due to the abundance of Z boson KK modes. 

The paper is organized as follows: In Section 2, we present the theoretical expression for the 
decay width of the LFV decay r — > iiV^v^ % — e, //, r, in the framework of the model III, with 
the inclusion of one and two NUEDs. Section 3 is devoted to discussion and our conclusions. 
In the appendix section, we give explicit expressions of the functions appearing in the general 
effective vertex for the interaction of off-shell Z-boson with a fermionic current. 

2 t — ► [i V{ V{ decay in the general two Higgs doublet 
model with the inclusion of non-universal extra di- 
mensions 

The LFV r — > /i Vi v { decay exists at least in the one loop level and, therefore, the physical 
quantities like the BR contains rich information about the model used and the free parameters 
existing. The extension of the Higgs sector in the SM makes the flavor violation (FV) possible 
with the help of the new Yukawa interactions coming from the new Higgs scalars. In the multi 
Higgs doublet models, the flavor changing neutral current (FCNC) at tree level, which induces 
the FV interactions, can appear and the general 2HDM, so called model III, is one of the 
candidate. The additional Yukawa interactions arising from new Higgs doublet are responsible 
for the non vanishing theoretical values of the physical quantities like the BR. The inclusion of 
extra dimensions may bring further contributions to these physical quantities. In the present 
work we assume that the first Higgs doublet and the gauge fields feel the extra dimensions, 
however, the leptons, the quarks and the second Higgs doublet, which contains new Higgs 
particles, live in the 4D brane. With the addition of single (double) extra dimension, the 
Yukawa interaction, including the LFV part in the model III, reads 

£>Y = vi(6)ijkL(<t>l\y(z)=o)E jR + C///, 7. 0-J /'-';/.' + h.C , (1) 

where i, j are family indices of leptons, L and R denote chiral projections L(R) = 1/2(1 =f 75), 
<f>i for i — 1, 2, are the two scalar doublets, Ul and Ejr are lepton doublets and singlets respec- 

2 



tively. Here vi(6)ij are 5(6)-dimensional dimensionful, ^ are dimensionless Yukawa couplings 
and vi(6)ij can rescaled to the ones in 4-dimension as vi(6)ij = V2ttR (2ttR) rj^, R is the 
compactification radius and y (z) is the coordinate represents the 5(6)'th dimension. The cou- 
pling which has complex entries in general, is responsible for the FCNC at tree level 1 . At 
this stage we take two Higgs doublets 4>i and (fi 2 as 

respecting that only the first one has a vacuum expectation value but not the second, namely. 



\ , ( V2X 
v + H° ) + { t X ° 



^2\H 1 + iH 2 



(2) 



<0i>=^(°) ;<0 2 >=o. (3) 

Therefore the CP even neutral Higgs particles Hi and H° do not mix in the tree level and they 
are obtained as the mass eigenstates h° and H° respectively. The neutral Higgs particle H 2 is 
the well known CP odd A°. In this case, the SM particles lie in the first doublet and the new 
particles in the second one. 

The LFV r — > [iViVi decay is induced by r — > fiZ* transition and Z* — > Vi Vi process 
(see Fig. [IJ. The r — > fiZ* transition, which needs the FCNC at tree level, is driven by the 
internal new neutral Higgs bosons h° and A , which are living in the 4D brane. Since the gauge 
fields are accessible to the extra dimensions, the internal Z boson has KK modes, which have 
additional contribution to the physical quantities related to the decay studied. The KK modes 
of gauge fields appear after the compactification of the external dimensions on the orbifold 
S 1 /Z 2 ( (S* 1 x S 1 )/Z 2 ) for a single (two) extra dimensions and, for two extra dimensions, the 
gauge fields can be expanded to the KK modes as: 

A,(x,y,z) = 7 ^^AAf Q \x) + 2^f j A^\x)cos{ny/R + Tz/R)\ , 



{27TR) d / 2 



n,r 



A*{x,y,z) = ^^ 2 ^Y,^ r \x)^Hny/R + rz/R)y (4) 
where % = 5, 6, d = 2, the indices n and r are positive integers including zero, but both are not 



zero at the same time. The KK modes of the gauge field Z have the masses \jvr? z + vn? n + m^, 
where m n = n/R and m r = r/R. Notice that the gauge boson mass matrix is diagonal since 
the Higgs field, which has the non-zero vacuum expectation value, is in the bulk (see jSU] for 
details). Here, we take the compactification radius R as the same for both new dimensions. In 
the case of a single extra dimension, one should set d — 1, take z — 0, and drop the summation 
over r in eq. (Jl|). 

1 Notice that, in the following, we replace £ E with ^ where "N" denotes the word "neutral". 
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Now, we present the general effective vertex for the interaction of off-shell Z-boson with a 
fermionic current 

T (REN) (r _> ^ = h lft + f 2 7^75 + h <V*" + h <V7s^ , (5) 

where k is the momentum transfer, k 2 = (p — p') 2 , p (p') is the four momentum vector of 
incoming (outgoing) lepton and the explicit expressions for the functions fx, f'2, fz and f'4 are 
given in the appendix section. The matrix element M of the r — > /x Vi Vi process is obtained 
with the internal Z boson connection between the r — > \x transition and the V{Vi pair. Since 
the gauge fields are accessible to the extra dimensions, the KK modes of the internal Z boson 
have additional contributions to the process (see Fig. Notice that the KK number need 
not to be conserved in the vertices where Z boson appears since the the extra dimensions in 
this case are so called NUED. Using the matrix element M the decay width T of the decay 
under consideration can be obtained in the r lepton rest frame with the help of the well known 
expression 

rfr ^iM Wp -|>)n^, (6 ) 

where p (pi, i=l,2,3) is four momentum vector of r lepton (/z lepton, incoming u, outgoing v). 

3 Discussion 

The r — > /j z/j is induced by the LFV r — > /i transition which depends on the various Yukawa 
couplings ^ 2 i = fi,r; j = e,fi,r and they need to be restricted by using the present 
and forthcoming experiments. In our work, we take only the r lepton as an internal lepton 
and, therefore, we choose the couplings £§ tTT and £§ T „ as non-zero. Here we expect that the 
couplings which contain at least one r index are dominant similar to the Cheng-Sher scenario 
|S2] • The upper limit of the coupling has been estimated as 30 GeV (see and references 
therein) by assuming that the new physics effects can not exceed experimental uncertainty 10~ 9 
in the measurement of the muon anomalous magnetic moment. In our numerical calculation 
we choose £§ tTfl = I GeV by respecting this upper limit. Since there is no restriction for the 
Yukawa coupling £jv rr , the numerical values we use are greater than £jf T „. 

This work is devoted to the a single and two NUED effects on the BR of the LFV processes 
r — > iiViVi y in the type III 2HDM. In the case of two extra dimensions, we observe that 
the contribution of KK modes enhances the BR considerably, due to the crowd of Z boson KK 

2 The dimensionfull Yukawa couplings are defined as — J £y y 
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modes. Notice that we use the numerical values, m z = 91 (GeV), m w = 80 (GeV), s w = v^0-23, 
G F = 1.6637 x 1(T 5 (GeV~ 2 ) } T T = 2.27 x 10~ 12 (GeV), in our numerical calculations. 

In Fig. |21 we present £§ TT dependence of the BR for £x >Tfl = 1 (GeV) and for different values 
of the compactification scale 1/R, in the case of a single extra dimension. Here solid (dashed, 
small dashed, dotted, dash-dotted, three dashed-spaced) line represents the case without extra 
dimensions (with extra dimensions for 1/R = 200, 400, 800, 1000, 2000 (GeV)). It is observed 
that the BR lies in the range 1.0 x 10~ 6 — 5.0 x 10~ 5 for the interval of the Yukawa coupling 
10 (GeV) < £jv TT < 50 (GeV). The addition of a single extra dimension enhances the BR 
almost twice for the small values of the compactification scale. However this enhancement is 
not more than the order of 1% for its large values. 

Fig. El represents the compactification scale 1/R dependence of the BR for £§ tT/Ji = 1 (GeV) 
and for different values of the coupling £,§ TT , in the case of a single extra dimension. Here solid- 
dashed-dotted straight lines (curved lines) represents the BR for £§ TT = 10 — 30 — 50 (GeV) 
without extra dimensions (with extra dimensions). This figure also shows that the contribution 
due to the extra dimensions is negligible especially for 1/R > 600 (GeV). 

Now we make the same analysis for two NUED. In this case the crowd of Z boson KK 
modes cause to enhance the BR of the decay analyzed and these effects can be observable. 
Notice that there is a possible divergence problem due to the abundance of KK modes in the 
summation of the KK mode contributions, however, the ratio m 2 + ( n 2 +r 2y R 2 appearing in the 
internal line converges to zero sharply with the increasing values of the integers n and r and 
the convergence of the KK sum is obtained for the region of the compactification scale that we 
study, 1/R > 200 GeV. 

In Fig. HI we present £§ >TT dependence of the BR for £^ rM = 1 (GeV) and for different values 
of the compactification scale 1/R in the case of two NUEDs. Here solid (dashed, small dashed, 
dotted, dash-dotted, three dashed-spaced) line represents the case without extra dimensions 
(with extra dimensions for 1/R = 200,400,800, 1000, 2000 (GeV)). This figure shows that the 
BR is considerably enhanced for the small values of the compactification scale 1/R. The BR 
is two orders (one order, three times, two times, 30 %) larger for 1/R = 200 (GeV) (1/-R = 
400 (GeV), 1/R = 800 (GeV), 1/R = 1000 (GeV), 1/R = 2000 (GeV)) compared to the one 
without NUED. Even for large values of the compactification scale near 1/R ~ 2000 (GeV) 
there is an enhancement in the BR. 

Fig. El is devoted to the compactification scale 1/R dependence of the BR for £j| T „ = 
1 (GeV) and for different values of the coupling £,§ TT , for two NUEDs. Here solid-dashed- 
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dotted straight lines (curved lines) represents the BR for Cn,tt = 10 — 30 — 50 (GeV) without 
extra dimensions (with extra dimensions). The enhancement in the BR is observed in this 
figure also. These contributions become negligible for the large values of the compactification 
scale, namely for 1/R > 2000 (GeV). 

At this stage we would like to summarize our results: 

• We predict the BR in the range 1.0 x 10~ 6 — 5.0 x 10~ 5 for the interval of the Yukawa cou- 
pling 10 (GeV) < £n,tt < 50 (GeV). The inclusion of a single NUED brings contributions 
which increase the BR twice of the one without the extra dimension, for the small values 
of the compactification scale. However this enhancement is not more than the order of 
1% for the large values of the compactification scale. 

• With the inclusion of two extra dimensions the BR is considerably enhanced for the small 
values of the compactification scale l/R, even two orders. In the case of large values of 
the compactification scale, there is still an enhancement in the BR. 

Therefore, the future theoretical and experimental investigations of the process r — > /i z/^z/j 
would ensure a valuable information about signals coming from the extra dimensions. 
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4 Appendix 



The explicit expressions for the functions fx, /a, and f± appearing in eq. (jSJ) read 

h = ra 2 a I dx ^— T \c v (m h +m h ) 

64 7T 2 cos 9 W Jo mf 2 — mf x \ 

(rself j^self 
(-mi r]f + m h (-l + x)7]J) ln—^- + (rrii r]f - m h (-l + x) r]J ) In-^f- 

+ {vriirjl + m h (-l + x) rtf) In l 'f - (77^77.+ + m h (-l + x)r)Y) In 2, f j 

K ' A* / 

+ c A (m h -m h ) 

(j^self j^self 
(-rrii r]~ + m h (-1 + x) rjf) In -^f- + (m 8 T]~ + m h (-1 + x) rjf) In 
/i 2 - fx 2 

+ [rriii]' + m h (-l + x) nf) In l, f + {-mini + mi 2 {-l + x)r]f) In 2, f \ 

a 2 ~ ix 2 ) J 

' dx / dylm 2 ( CAV f-c v r)Y)(-^- + — -) 



64 vr 2 cos 9 W Jo Jo » { lK A H y L\% Lft 

(l-x-y)mi [c A (m h - m h ) rfc (-— - + c v (m h + m h ) r)f (-^ + -^)) 

(Y Y J^ver jjver 

- 2 + (k 2 x y + m h m h (-1 + x + yf) (-— + — : ) - Zn -f- ^f- 
-^feo -^a A* A* 

(r]f{xm h +ym h )+m i r]l r]f(xm h +ym k )-m i r]l\ 
(m h + m h ) (1 - x - y)\ — — + 



O T ver o r ver 

/ ^AO h ( > A AO 



1 j ver j ver "\ 

+ 2^ /i 2 /i 2 J ' 

•^ 2 = RA 2 a I dx ^— 2\ C v{ m h- m h) 

64 tc 2 cos W Jo mf 2 — m 2 h \ 



j^self j^self 

(mi 7]~ +m h (-l + x) r]f ) In l 'f + (-m, + m h (-l + x) rjf) In 2 'f 

j^self j^self \ 

+ {-mini + m h (-l + x) rjf) In-^f - + (rrii n~ + m fa (-l + x)r]f) ln-^f- 
+ c A (m i2 + m h ) 

Cj^self j^self 
(mi r& + m ;i (-l + x) r/f) in ^ - (m 8 r/+ + m/ 2 (-l + x) iff) In 2 ' A ° 
/i 2 - /x 2 

+ (-mi nf + m h (-l + ar) r/f ) -^f - + (m 8 r/+ - m h (-l + x) nf) 

/i fi 

+ 7T- — 7, — i dx [ dyl m 2 (cy rif — c A vY) ( 1 ) 

64ti 2 cos6 w Jo Jo y \ tK h U 3 ^ L v ^> 
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- rriiil-x-y) (c v (m h - m h ) ^ + c A (m h + m h ) rjA (-^ - 7^7) 

(■y I jjver j^ver 1 

- 2 + (k 2 x y - m h m l2 (-1 + x + yf){-— + -—) -In^^f- 

_ w _ s ( rjj_ (xm h -ym h ) + | rtf (xm h - ym h ) - m i r)t\ 

\ m h m h) \ L x V) \ o rver ' 



o r ver o r ver 



i r ver T ver 

_ -rj v In A ° h ° h ° A ° 



2' 1 /i 2 /i 2 J ' 

= -i TA 2 9 a I dx I d y\\( 1 - x -y)( c vvY + c A r)f)(xm h +ym l2 
64 7T Z cos Uw Jo Jo V 



64 7T 2 COS 6*^ 

+ m i (c j4 (a;-y)^r + c v^ + (^ + y)) 



r ver 



+ (^(1 - x - y) (c v rjY + c A r)f) {xm h +ym h ) - rrii (c A (x - y)^ + c v rfc (x + y))j 
(1 x J vt{xm h +ym l2 ) , 1 1 x m lV r , 1 

v "/I 9 V Tver Tver J o V Tver Tver I I | ' 

h = ~ i 6 4 7r 2^ os / rfa; / ^{ U 1 ~ x ~y) {- ( c vvt + cavY) (xm h -ym h )) 



T ver 
L A° 



rrii (c A (x - y) rfc + c y ^ (x + y)) 



r ver 



+ ^(l-z-y)(-(c y ?7f + c A 77f)(:cmz 1 -ym h )) + rrii (c A (x - y) r]+ + c v r] { (x + y))^j 



where 



r ver T ver J o \ T ver r ver 

^A°h° ^h»A» A ^A°h» ^h°A° 



L i% = m l° (1 - x) + {m 2 - m 2 x (1 -x))x, 

T sel f T sel f i \ 

L^o = L 1 J {m h o -> m A o) , 

r self T self / \ 

L 2,h° = L i,ho{m h ->m, 2 ), 

rself T self i \ 

L 2,A° = L l,Ao{m h ->7Tli 2 ), 

^M» r = (! - x - y) + m\ (x + y)-k 2 xy, 

LhPAQ — m A0 x + m i (1 - x - y) + (m 2 h o - k 2 x) y, 

L v A e S = L v ^ r (m h o -> m A o) , 

L^ e o\o = Iw AO (m h o ^m A o), (8) 



and 



Hi — <,N,hi<,N,il2 Z,N,ih^N,l 2 i ' 
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A _ pE pE* _ pE* pE 
'li — $N,h&N,U2 ^N,ih^N,l 2 i ' 

+ _ pE* pE* , pE pE 
-B* tE* fE tf _E 



The parameters cy and C4 are ca = — \ and cy = \ — sin 2 9w in eq. © the flavor changing 



couplings ^ j represent the effective interaction between the internal lepton i, (i = e,fi,r) 



and outgoing (incoming) j = li (j = l 2 ) one. Here we take only the r lepton in the internal 
line and we neglect all the Yukawa couplings except £§ TT and £^ )T „ in the loop contributions 
(see Discussion section). The Yukawa couplings £|| ^ are complex in general, however, in the 
present work, we take them real. 
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Figure 1: One loop diagrams contribute to r — > /i z7 i z/j, i = e,/i,r decay due to the neutral 
Higgs bosons /io and A in the model III version of 2HDM. Solid lines represent leptons and 
neutrinos, curly (dashed) lines represent the virtual Z boson and its KK modes in 6 dimensions 
(h and A fields). 
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Figure 2: £jv rT dependence of the BR for £§ T „ = 1 (GeV) and for different values of the 
compactification scale 1/R, in the case of a single extra dimension. Here solid (dashed, small 
dashed, dotted, dash-dotted, three dashed-spaced) line represents the BR without extra dimen- 
sions (with extra dimensions for 1/R = 200, 400, 800, 1000, 2000 (GeV)). 
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Figure 3: The compactiflcation scale l/R dependence of the BR for £§ Tfl = 1 (GeV) and for 
different values of the coupling £jy TT , in the case of a single extra dimension. Here solid-dashed- 
dotted, straight lines (curved lines) represents the BR for ^§ TT = 10 — 30 — 50 (GeV) without 
extra dimensions (with extra dimensions). 
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Figure 4: £§ TT dependence of the BR for £§ Tfl = 1 (GeV) and for different values of the com- 
pactification scale 1/R, in the case of two extra dimensions. Here solid (dashed, small dashed, 
dotted, dash-dotted, three dashed-spaced) line represents the BR without extra dimensions 
(with extra dimensions for 1/R = 200, 400, 800, 1000, 2000 GeV) 
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Figure 5: The compactiflcation scale l/R dependence of the BR for £§ Tfi = 1 (GeV) and for 
different values of the coupling £jy TT , in the case of two extra dimensions. Here solid-dashed- 
dotted, straight lines (curved lines) represents the BR for ^§ TT = 10 — 30 — 50 (GeV) without 
extra dimensions (with extra dimensions). 
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